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ABSTRACT 

Using the 3.6/4.5/5.8/8.0 images witli 1.2" pixel resolution from the Spitzer/GLIMPSE 
survey, we investigated 23 masing and 18 very red objects that were not identified in the 2MASS 
survey. Counterparts for all selected objects were found in the GLIMPSE images. Color indices 
in these IR bands suggest the presence of a high-extinction layer of more than a few tenths of 
a solar mass in front of the central star. Furthermore, radio observations in the SiO and H2O 
maser lines found characteristic maser-line spectra of the embedded objects, e.g., the SiO J — 1-0 
line intensity in the v = 2 state stronger than that of the v = \ state, or very widespread H2O 
maser emission spectra. This indicates that these objects are actually enshrouded by very thick 
circumstellar matter, some of which cannot be ascribed to the AGB wind of the central star. 
Individually interesting objects are discussed, including two newly found water fountains and an 
SiO source with nebulosity. 

Subject headings: infrared: stars — stars: AGB and post-AGB — stars: mass loss — radio lines: stars 
1. Introduction to the realization that the maser emissions are 



J -a, ma nifestation of the mass loss process (|Kwok 

The discovery of stellar maser sources (|Wilson fc Barret fj gjgD ^ ^nd OH/IR stars are hkely to be stars in 

\mm and their association with very red infrared t^e asy mptotic giant branch (AGB) phase of evo- 

objects has led to a new class of objects called i^^io^ (lHabing||l996l : ISevenster et all I2OO2I) . The 

OH/IR stars. From the OH profiles, we came ^^^ss loss process also creates a thick circumstel- 



lar dust envelope which often obscures the central 
star. Although O H/IR stars are often brig ht mid- 



infrared objects (jJimenez-Esteban et al. 



many of them have no optical or near infrared 
counterparts. In the absence of a spectral clas- 
sification based on photospheric spectroscopy, 
the nature of these stars is classified based on 
the circumstellar spectra in the mid-infrared 
(|Kwok. Volk fc Bidelma^ll997f ). 

It should be noted that the mid-infrared col- 
ors alone are not enough to determine the nature 
of maser sources. For example, the IRAS col- 
ors of mass-losing evolved stars overlap with col- 
ors of young stellar ob jects embedded in molec- 
ular clouds (jvan der Veen fc Habing . .1988,1 . and 



there is no well-defined boundaries in the color- 
color dia gram between young and evolved stars 
(jWeintraub , 1990). The maser characteristics by 



themselves are not enough to dis tinguish the two 
classes of objects ( Caswell 1999[) . The situation 
for SiO maser sources is also complex. For exam- 
ple, there ar e three SiO ma s er sources in molec- 
ular clouds (jHasegawa 1985 : Morita et al. Ill992 ) 
that are likely to be young stellar objects. There 
are also peculiar cases, such as the SiO maser 
source IRAS 19312-1-1950 is an evolved star in a 



dark c loud (jNakashima et al. 2004: .Deguchi et al 



2004d), and V838 Mon is an M supergiant with 



SiO emission which emerged after nova eruption 
(Deguchi et al. 2005b; Tylenda & Sokcr 2006). 

In some cases, the imaging and identification 
of the central star in NIR bands can greatly fa- 
cilitate the determination of the nature of thes e 
objects (jLewis et al. I l2004t iDeguchi et al]l2005a[ ). 
Young, low-mass stars are usually born in a star 
cluster, but evolved stars, especially AGB/post- 
AGB stars with the ages of Giga years, are ob- 
served as a single object (see § 3 of IPeguchi et al.l 
However, for the case of supergiants, they 
are frequently a m ember of stellar associations 



( Humphrevs 1970f) . suggesting they are on their 



way towards leaving the originating star clusters. 
Deep infrared imaging can help find the surround- 
ing low-mass stars. 

The recently performed Galactic Legacy In- 
frared Mid-Plane Survey Extraordinaire (GLIMPSE) 
of the Spitzer Space Telescope (SST) has provided 
unprecedented deep near-infrared images of the 
Galactic plane. Because of the high angular res- 
olution 1.2" per pixel) of the GLIMPSE observa- 



tions, we have the opportunity to obtain an accu- 
rate identification of the infrared/maser sources. 
In this paper, we report the identification of 41 
infrared sources from the GLIMPSE data, as well 
as SiO or II2O maser observations of these sources 
with the Nobeyama 45-m telescope. 

2. Observations 
2.1. Source selection 

We selected 41 sources from the Midcourse 
Space Experiment (MSX) Galactic plane survey 
that are found to have no (or dubious) near- 
i nfrared counterparts i n the 2MASS sky survey 
2000f) . Half o f these (22 ob 



( Skrutskie et al. 



j ects) a re OH/IR objects listed bvlSevenster et al 



( 2001 ). Ite Lintel-Hekkert et al. I (|l99lD . lEder et al 



19881), etc. The positional accuracy of the MSX 



posit ions is about a few arcseconds (jEgan et al 



19991) . comparable to the positional accuracies of 
OH masers measured with the Very L a rge A rray 
(VLA) [see figure 9 of iDeguchi et al. I (|2002l ) for 
comparison]. Most of these objects have reason- 
ably red mid-IR colors with Cce [= log{FE/ Fc)] 
between and 0.5, where Fc and Fe are the 
MSX band C (12 /im) and band E (21 ^m) flux 
densities, respectively. 

In addition to the MSX sources, we have 
included a number of medium-bright and red 
IRAS sources without 2MASS counterparts for 
which no maser line search was made. These 
are IRAS 18030-1858, 18151-1713, 18241-1010, 
18298-0904, 18424-0509, 18444-0359, 19011-1-0638, 
19025-1-0702, 19087-1-1006, 19112+1220, 191144-0920, 
19180-1-1230, 19254+1724, 19266+1815, 19275+1720, 
G014.2758-00.1602, G027.8433-00.1060. The 
MSX colors of these objects are not much dif- 
ferent from those of the above OH/IR sources 
without NIR counterparts. The entire sample we 
chose for identification in the Spitzer/GLIMPSE 
survey is listed in table 1, giving the source 
name, MSX(6C) name, right ascension and decli- 
nation (J2000) of the MSX source, flux density 
in the MSX C (12^m) band, the MSX colors 
[Cac = logiFc/FA), and C'ce = logiFE/Fc)], 
status of detections in SiO, H2O, and OH masers, 
and a conventional OH object name. 
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2.2. Spitzer/GLIMPSE Identification 

The GLIMPSE survey was made with SST be- 
tween March and November, 2 004 using the In- 



frare d Array Camera (IRAC) (jlndebetouw et al 
2004 ). The survey covered the area of\l\ = 10-65° 
over latitudes |5| < 1° toward the inner Galactic 
disk excluding the Galactic center. Simultaneous 
5.2' X 5.2' images at 3.6, 4.5, 5.8 and 8 /im were 
obtained at a spatial resolution of ^ 1.2" x 1.2" 
per pixel. 

The color-composite images (from the 3.6, 4.5 
and 5.8 fim band images) toward the 41 selected 
objects are shown in Figure 1. The counterparts 
for all sources in our sample were found; they are 
all very red and near the center of the images. 
Because all of these objects are bright in 8 fim, 
the identifications in other bands were made quite 
easily. The color images were created with special 
software developed at the Space Astronomy Labo- 
ratory of the University of Calgary. The flux den- 
sities were derived from the GLIMPSE archives; 
the aperture size used for photometry was typi- 
cally a few to 12 arcseconds, depending on the 
image size of the star. The estimated error of 
the photometry is dependent on the background 
noise levels, and are typically < 2% for the ob- 
jects brighter than 1, 3, 7, 16 mJy at the 3.6, 4.5 
and 5.8 /^m bands respectively. The flux densities 
are given in columns 4-7 in table 2. The color in- 
dices, [3.6] — [5.8] and [3.6] — [12], are also listed 
in the 8th and 9th columns in table 2. Here, [3.6], 
[5.8], and [12] are the magnitudes in the GLIMPSE 
3.6 and and 5.8 /im, and MSX 12 fim bands, re- 
spectively, and the flux densities at 0th magnitude 
are 280.9, 115.6, and 26.4 Jy in these bands, re- 
spectively. The interstellar extinction toward each 
object is given in the 10th column of table 2 and 
it is estimated by the formula based on the full- 
sky 100 fim composite map of COBE/DIRBE and 
IRAS maps (Schlegel et al. 1998)0 

After getting accurate GLIMPSE positions 
(the 2nd and 3rd columns in table 2), we again 
checked the 2MASS images and found the red 
counterparts corresponding to the positions of 
4 sources; 18286-0959 ( J18312292-0957194), 
18298-0904 (J18323700-0902301), 19275-1-1902 
( J19294645-hl709013), and 19266-1-1815 ( J19285303 
-1-1821228). These counterparts are very faint 



in the K-band {K — 12.6-14.7), making pre- 
vious identifications uncertain. Their positions 
coincide wit h those of the G LIM PSE o bjects 
within 0.5". Jimenez-Esteban et al. I (1200(1 listed 
NIR counterparts for 5 objects in the present 
sample: 18182-1504, 18460-0254, 18488-0107, 
19087-fl006, and 19254-^1631. Because of the 
faintness of the objects {K > 13), their identifica- 
tions were again doubtful. We have checked the 
GLIMPSE positions of these objects and found 
that 18182-1504, 18460-0254, and 19254-1-1631 
were correctly identified, but 18488—0107 and 
19087-H1006 were misidentified; the GLIMPSE ob- 
jects are located 5"NE and 6"S ofthe 2MASS stars 
(K = 13.05 and 14.10) given bv lJimenez-Esteban et al 
(l2006l) . 

In addition to the above selected objects 
without 2MASS counterparts, we checked the 
GLIMPSE images of about 200 SiO maser sources 
within the area of I = 10 - 60°, and |6| < 1° 
where the near-infrared identification s are a lready 
available (|Deguchi et al.1ll998ll200lll2002l ). The 
GLIMPSE images mostly exhibit a single bright 
object at the center of the field, except for IRAS 
19027-1-0517, which shows accompanying nebu- 
losity (Fig. 6). Because this object has a NIR 
counterpart in 2MASS images and seems to be 
intrinsically different from the objects in our sam- 
ple, we gave a discussion of this object separately 
in Appendix A. 

We wish to make a special reference to IRAS 
18450-0148 (W43A), which is known as a wa- 
ter maser fount ain with coUimated outflows 
(|lmai et al.ll2002h . The detection of SiO maser 



emissi on in this object ([Nakashima fc Deguchi 



^ available a.t \http://ned'ww'w.ipac.caltech.edu/forms/calculator.html\ 



2003br ) suggests that this is an evolved star. The 
near-infrared counterpart of this interesting object 
is identified here for the first time. 

2.3. SiO and H2O maser observations at 
Nobeyama 

Because a majority of objects in our sample 
have not been searched for SiO or II2O maser 
emissions before, we observed the objects in SiO 
J = 1-0 V = 1 and 2 and H2O 6i6-523 maser lines 
with the 45-m telescope at Nobeyama during April 
10-23, 2006. A cooled HEMT receiver (H40) and 
an accousto-optical spectrometer array were used 
for the 43 GHz observations. The receiver system 
temperature was about 200 K and the effective 
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velocity resolution is 0.3 km s ^. The half-power 
beam width (HPBW) at 43 GHz was ~40" and 
the conversion factor of the antenna temperature 
to the flux density was ^2.9 Jy K~^. In addition 
to the SiO maser observations, we made a 22.235 
GHz H2O maser observation on April 20 and 23, 
when weather was unfavorable for 43 GHz obser- 
vations. We used a cooled HEMT receiver (H22) 
and the same accousto-optical spectrometer array 
(with an effective velocity resolution of 0.6 km 
s~^). The conversion factor of the antenna tem- 
perature to the flux density was ^^2.8 Jy K^^ at 
22 GHz. Because the H40 and H22 receivers have 
a 2 GHz band width, we configured the AOS-H 
spectrometer array to have the spectra of the SiO 
J = 1-0 w = (43.423 GHz), and w = 3 (42.519 
GHz), 29SiO J = 1-0, v=0 (42.880 GHz), and 
H53a (42.952 GHz) lines simultaneously in the 
H40 receiver, and NH3 Jk — li-li, 22-22, and 
33-33 lines (23.694, 23.722, 23.870 GHz, respec- 
tively) simultaneously in the H22 receiver. The 
rms noise temperatures for these additional lines 
are similar to the noise levels of the SiO and H2O 
lines and therefore are not repeated in tables 3 and 
4. No detections in these additional lines were ob- 
tained except for 18182—1504. Further details of 
observations using the N RO 45-m telescope are 



describ ed elsewhere (e.g., iNakashima fc Deguchi 
2003ah . 



Among the 31 objects observed, 8 were detected 
in SiO and their spectra are shown in Fig. 2. In all 
cases, the v = 2 maser line is stronger than the v = 
1 maser line. Water masers in the 6i6-523 tran- 
sition were detected in two sources (18286—0959 
and 18460—0151) and their spectra are shown in 
Fig. 3. It is notable that both of the detected 
objects exhibit a wide velocity range in the H2O 
maser emission spectra, which res emble th e H2O 
spectra of "water fountains" (Mor ris et al ][2003). 

Because some of our sample objects might be 
associated with young stellar objects in molecular 
clouds, we also made simultaneous observations of 
the HCO+ J = 1-0 hne at 89.189 GHz and the 
SiO J = 2-1 V = 1 line. The non-detection of the 
HCO+ line effectively excludes the possibility of 
these water maser sources being associated with 
young stellar objects. 



3. Discussion 

3.1. Tvifo-color diagram 

The selected objects have extremely thick cir- 
cumstellar envelopes. In addition, they suffer from 
a large interstellar reddening in the Galactic plane 
even at NIR to MIR wavelengths. Therefore, it is 
hard to separate the circumstellar from the inter- 
stellar reddening for these objects because of the 
uncertainty of the interstellar extinction correc- 
tionll We discuss here quantities of the sampled 
objects by introducing the color indices, [3.6] — [5.6] 
and [3.6] — [12], which vary relatively mildly with 
the large extinction. We apply a wavelength- 
dependent extinction derived from the GLIMPSE 
survey which was obtained fro m the color excess 
of stars at I = 42° and 284° ( Indebetouw et al 



2OO5I ). We compare physical quantities of the sam- 
ple objects with those of the proto-typical OH/IR 
star with SiO masers, OH 127.8+0.0, which is be- 
lieved to be at the final stage of evolution on the 
AGB. Physical par ameters of this sta r are rela- 
tively well known (jSuh fc Kim I l2002l ) , with the 



optical depth of the envelope being about 10-15 
at 9.8 Aim. This bright OH/IR star [Fc = 184 
Jy) has [3.6] - [5.8] = 3.46 and [3.6] - [12] = 5.0. 
Because of the position {I ~ 128°), the interstellar 
extinction at the NIR wavelengths is small for this 
object [Ak = 0.64). If this object were behind a 
dust cloud with a heavier extinction of Ak > 7.6, 
it would be undetectable on the 2MASS image 
[K > 13.7). The sphe rically symmet ric model for 
OH 127.8+0.0 (Ke mper et al. II2OO2I ) gave a mass 
loss rate of 7 x 10~^ Mq yr~^. The same model 
gives a total mass of the envelope of 0.2 M© (if 
integrated to the outer r adius of the dust she ll, 
1.3 X 10^^ cm, as used bv lKemper et al. I (|2002l )). 

Figure 4 shows a two-color diagram, [3.6] — [12] 
versus [3.6] — [5.8], for the sampled objects, where 
OH 128.7+0.0 occupies the position indicated by 
the square in this diagram. The interstellar red- 
dening moves the position of the star to the upper 
right with an inclination of about unity because 
the extinct ions at 5.8 an d 12 fj,m are comparable 
(Table 1 of iMathis Ill990l ). Selected objects in this 
diagram distribute from the lower-left to the up- 



^ As noted in Section 2.2, we used a formula given by 
ISchlegel et "aD (1998), which gives almost the maximum 
value of interstellar extinction due to thick dust clouds to- 
ward each source. 
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per right with a steeper slope (~ 1.71 ± 0.19) than 
the slope of the interstellar reddening line (< 1). 
They appear in a very wide range of color indices, 
over 3 in [3.6] — [5.8] (corresponding to Ak = 20), 
which is too large to be considered an effect of in- 
terstellar reddening. These facts suggest that the 
extremely red colors of these objects cannot be ex- 
plained solely by the interstellar reddening applied 
to the OH 127.8+0.0 type star, but they are more 
or less intrinsic to these objects. 

The steepness of the distribution seems to be 
understandable by the model that the 12 /im flux 
is a re-emission of absorbed NIR light by dust 
grains near the central star. If MIR re-emission 
does not occur (as interstellar reddening) , the star 
should move on the line indicated by the dotted 
arrow. SiO detections (filled circles) also appear 
in a wide range of color indices in figure 4, but 
the SiO sources distribute with a steeper inclina- 
tion (~ 1.98 ± 0.35) than the no-SiO subsample 
1.62 ± 0.21), indicating that the re-emission ef- 
fect at 12 ^m is stronger in the SiO maser sources 
than the non-SiO emitting objects. The SiO maser 
is an indicator of active mass loss near the central 
star (^ 10^'* cm), suggesting the presence of thick 
dust re-emitting the stellar radiation more at MIR 
wavelengths than the objects without SiO masers. 

The above findings strongly suggest that these 
objects have an excessively large optical depth of 
circumstellar dust, which is much larger than that 
of the circumstellar dust of OH 127.8+0.0. The 
order of magnitudes of the excess material can be 
estimated from [3.6] — [12]. We use the relation 
between [3.6] — [12] and the rn ass loss rate calcu- 
lated bv iKemper et al. I (2002). Applying the in- 
terstellar reddening corrections to [3.6] — [12], we 
computed the mass loss rate of the envelope from 
figure 6of lKemper et al. I(|2002h . and obtained the 
excess factor (/c) of the circumstellar matter (rela- 
tive to that of the OH 127.8+0.0 envelope), which 
is listed in the last column of table 2. They exceed 
the optical depth of the OH 127.8+0.0 envelope by 
a factor of a few up to 80. Uncertainty of the inter- 
stellar extinction correction of about 30 % makes 
a shift of [3.6] — [12] up to 0.7 at most, making 
the uncertainty of the excess factor to be within a 
factor of 2 even for the worst case. 

The total envelope mass for each sampled ob- 
ject must exceed by a similar factor, i.e. these ob- 
jects might have envelope masses of 0.4-20 Mq if 



simply integrated until the outer radius of ~ 10^'' 
cm. The total mass of the envelope is propor- 
tional to the outer-boundary radius for the as- 
sumed density distribution of a constant-mass-loss 
model (p oc r~^), whereas the optical depth is 
rather sensitive to the inner-bounday radius of the 
dust envelope. Because the color index, [3.6] — [12], 
which was used here for evaluation of the excess 
factor, does not give useful restrictions on the 
outer boundary radius in the Kemper et al's mod- 
eling (it is rather restricted by IRAS 60 or 100 
lim flux density), the total mass of the thick mat- 
ter estimated here has large ambiguity. Because, 
without precise modeling, it is uncertain as to how 
far from the central star the extra material in the 
envelope is located, we only use the excess factor 
fc in this paper, which characterizes the envelope 
of these objects. 

The total envelope mass for each sampled ob- 
ject can be reduced, in some degree, by introduc- 
ing a non-spherical distribution of the dust enve- 
lope, e.g. a torus or disk structure. However, even 
a wide opening angle of 120° for cavity cones can 
reduce the mass by 50% compared with that of the 
spherical distribution with the same radial density 
profile. As discussed in the later sections, high- 
velocity components of water maser emission may 
suggest the presence of the accretion disk which 
creates collimated jets. In such a case, the scat- 
tered light through the pole of the torus should be 
observable, although it is hard to detect the scat- 
tered light for these objects because of the large 
interstellar extinction in the K band. 

3.2. Characteristic of SiO maser emission 

It is striking that only the SiO J = 1-0 v = 2 
line was detected in 6 out of 8 SiO detected 
sources. Even in the remaining two sources the 
intensity of the J = 1-0 v = 2 line was consid- 
erably stronger than that of the J — 1-0 v = 1 
line. The trend, i.e. the increase of the v = 2/1 
line intensity ratio with the infrared color [C12 = 
log{F25/Fi2), where F12 and F25 are IRAS 12 
and 25 ;[xm flux density], was clear l y first d emon- 
strated by iNakashima fc Deguchi (2003b). For 
the present sample of SiO detected sources, the 
MSX color, CcE = log{F2i/ F12), is between 0.10 
and 0.34, which fits well with the above trend. 
All of the SiO detected sources were previously 
observed by the OH 1612 MHz line and were all 
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detected except for 18241-1010. 

Figure 4 indicates that, for the sampled objects, 
SiO maser detection rates do not seem to correlate 
with the colors, [3.5] - [5.8], or [3.6] -[12]. The flat 
detection rate seems to suggest that these stars are 
still in a mass-losing stage, at the final transient 
stage of the AGB to planetary-nebula phase. They 
are not at a later stage of the post-AGB phase 
when SiO masers should disappear. 

3.3. New H2O sources with wide spread 
emission 

Two H2O maser sources, 18286—0959 and 
18460—0151 (Fig. 3) have a very wide velocity 
range (^200 km s^^) in their maser emission spec- 
tra. Although the rich H2O maser emission spec- 
trum in IRAS 18286-0959 (Fig. 3 left) resembles 
the water maser spectra a ssociated with compact 
H II regions (for example, Ku rtz fc Hofne r 2002), 
the nondetections of HCO+, H53q!, and NH3 lines 
toward these sources (the present work: table 3) 
suggest that they are not associated with dense 
molecular clouds (or compact H 11 regions), and 
the rich H2O emission could be the result of a very 
irregular velocity field in the envelope. The fact 
that these sources are point-like with no associated 
nebulosity in their GLIMPSE images suggests that 
they are stellar objects. 

Similar wide spread water maser spectra were 
found in the class of "water fountain" sources. 
There are four known water fountain sources 
[16342-3814, OH 12.8-f 0.9 (181 39-1816), W43A 
(18450-0148), 19134-H2131 dLjk kcl fc Morris 



1988t iBoboltz et al. |[20o5 llmai et al. 2004. 200 



with a few more p r obabl e cases ([Gomez et al 



19941: 

jects exhibit 



Deacon et al. I l200lh . 



All of these ob- 
maser components in 



the H2O 

a velocity range that exceed the OH velocity 
range. The H2O maser components are spatially 
more extended t h an the OH maser co mponents 
(jlmai et al.l 120041: iBoboltz et al. I [20051 ). Among 



these "water fountain sources", 16342—3814 had 
the highest velocity separation of ^ 258 km s^^. 
Our object 18460—0151 definitely belongs to the 
class of water fountain sources from the similarity 
of the emission characteristic, and it has a record- 
high separation of 292 km s"^ between emission 
components. 

Both new objects, 18286-0959 and 18460-0151, 



are relatively bright IRAS sources (25 and 20 Jy 
at 12 iim respectively); the former has a 2MASS 
counterpart, J18312292- 0957194, with K magni- 
tude of 12.67 and - if = 0.89. They have a ris- 
ing mid-IR spectrum toward longer wavelengths, 
indicating the presence of a thick dust envelope. 
From the radial velocities, we estimate kinematic 
distances of 18286-0959 and 18460-0151 as 3.1 
and 6.8 kpc, respectively, and luminosities as 
8.7 X lO'^ and 4.2 x 10^ Lq for these distances. 

If these objects have a thick dust torus (as in 
the unified model of the Type I and II active galac- 
tic nuclei), o ptically observabl e objects (as IRAS 
16342-3814; ISahai et al.lll999l ) must be seen from 
the polar-axis direction of the dust torus and un- 
seen from equatorial directions. If these high- 
velocity features are a part of the polar jet, which 
is created by the accretion disk, the higher velocity 
objects should have bluer colors. This hypothesis 
seems applicable for 18286—0959 (because of the 
identified NIR counterpart with K — 12.67), but 
not for 18460-0151 (with no NIR counterpart). 
Regarding 18460—0151, which has the highest ve- 
locity separation among the water fountains, the 
large distance (~ 6.8 kpc) and the strong interstel- 
lar extinction hide this object behind molecular 
clouds. Applying the large interstellar extinction 
correction of Ak = 8.2, we infer the extinction 
corrected K magnitude to be fainter than 5.8. It 
is slightly fainter than the extinction-corrected 3.6 
/im magnitude ([3.6] ~ 5.5), suggesting that the 
scattered light is not seen from this object. 

3.4. Envelopes of extremely red OH/IR 
stars 

The central stars of extremely red OH/IR ob- 
jects with low expansion velocity {Vexp < 15 
km s~^) have been considered to be relatively 
low-initial-mass post-AGB stars compared with 
those of normal-color sources (jSevenster et al.l 
2OO2I). In fact, t h e red ("RI") group sources 



of Sevenster et al. ( 20021 ) involved two objects 
in the present sample: IRAS 18135—1456 and 
18596-1-0315, which have the expansion veloci- 
ties, 14.8 and 13.6 km s""'^, respectively. Regard- 
ing expansion velocites derived from OH peaks, 
two water fountains, 18450-0148 (W43A) and 
18460—0151, are low expansion-velocity (lower- 
mass) sources (Vexp = 6 and 11 km s~^). 

Figure 5 shows a plot of the excess factor of the 
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envelope versus expansion velocity for OH 1612 
MHz doubly-peaked objects in the present sam- 
ple. Note that the excess factor was derived from 
[3.6] - [12] (relative to OH 127.8+0.0; see section 
3.1). It indicates that the excess factor does not 
correlate strongly with the expansion velocity of 
the envelope, and hence the initial mass of the 
central star. 

Figure 5 also indicates that objects with small 
expansion velocity {Vexp < 15 km s^^) seem to 
have lower SiO detection rates than the larger 
expansion- velocity objects {Vexp > 15 km s~^). 
This indicates that the lower-mass stars with small 
expansion velocities are at the phase unfavorable 
for making strong SiO masers, e.g. dissociat- 
ing SiO molecules by hot central star radiation, 
and/or terminating mass loss fr om the central 
star. This is consistent with the ISevenster et al 



( 20021 ) 's finding that the red objects with small 
expansion velocites are stars in the early post- 
AGB phases with nearly zero mass loss; the ex- 
treme high- velocity flow develops at this phase and 
SiO masers gradually disappear. We infer that the 
red stars with large expansion velocites are AGB 
stars still in an active mass-losing stage because of 
the presence of SiO masers. However both groups 
of stars, which are investigated in this paper, are 
dressed by excessively large amounts of circumstel- 
lar matter irrespective of the mass of the central 
star. 

The radii of the OH emitting regions were found 
to be 3 - 16 xlO^^ cm for some of the selected 
objects (OH 21.5-h0.5, 30.1-0.7, and 32.0-0.5) 
with the OH phase-lag / angular-size measurements 
(Herman et al. 19851). Because the HoO mas- 



ing region of water fountains i s extended more 
than the OH emission region ( Imai et al. 20021 



for W43A), the thick material must be extended 
out of the OH masing region. There is evidence 
that H2O maser outflow strikes the dense mate- 
rial which is located outside of the OH masing 
region. Therefore, it is likely that the excess ma- 
terial found in section 3.1 is extended to consid- 
erably outer parts of the envelope, i.e. outside of 
the OH masing region. The H2 number density 
of the OH masing regior i is known to be ^ 10'^ - 
10^ cm"3 (for example, iNetzer fc Knapp |[l987h . 
and the density of the H2O masing region must 
be, by several orders of magnitude, higher (~ 10^ 
cm~'^). Thus, we cannot deny the possibility that 



the thick layer is extended at the radius of more 
than 3 x 10"'^^ cm, though such a layer could be 
clumpy and dense, depending on the model. En- 
velope masses of some OH/IR stars [e.g., IRAS 
18450-0148 (W43A)], exceed 4 Mq even for the 
outer radius of 3 x 10^^ cm. Though the envelope 
mass can vary by a factor of a few, depending on 
the uncertainty of interstellar extinction and non- 
spherical distribution of the thick layer as noted 
before, it is difficult to lower the envelope mass 
less than 1 AIq for several of the thickest objects 
of the sample. 

A number of numerical modelings of thick cir- 
cumstellar envelopes in the outg oing AGB phase 



van Hoof et al. 



have been made (for example, , 

1997t Lorenz-Martins fc de Arauio Il997t ). iDavid fc Papoular 
(119921 ) and lSuh & Jones I (|l997h examined super- 



wind models to explain far IR spectra of OH/IR 
objects accompanying deep silicate absorptions. 
One of their models, which has a dense shell be- 
tween radii 2 x 10^^ and 4 x 10^^ cm in addition 



to a normal continuous flow (jSuh fc Jones Ill997f ). 



successfully fits the computed to the observed 
spectral energy distribution of IRAS 18257-1000 
(OH 21.5-1-0.5), which has a relatively mild excess 
factor (1.4) in our fist (table 2). The large IRAS 
60 /im flux densities of most of the other objects 
(though some of them might be contaminated 
from nearby clouds) entail for the much denser 
outer shell in these models, which inevitably re- 
sults in a large total envelope mas s of more than a 
few Mq (jPavid fc Papoular 111 992j ): the superwind 
models artificially restrict the outer radius of the 
thick material to keep the enclosed mass smaller. 
Although the total mass of the outer envelope 
could in principle be observationally constrained 
by CO J = 1-0 intensities or IRAS 60/100 ^im 
flux densities, it is hard to estimate for these ob- 
jects because of contamination by radiation from 
surrounding clouds. 

Though all of these modelings of very red 
OH/IR objects assume a superwind with a mass 
loss rate of 10~^ Mq yr~^ in a duration of 
about a few hundred years, the origin of the thick 
material is not necessarily restricted to the su- 
perwinds which occur at the end of the AGB 
phase. Because no correlation appears between 
the excess factor and the expansion velocity, 
the thick material does not seem to be related 
with the central star masses. The excess mate- 
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rial at the outer envelope of these objects can 
be a source of a shock front of the extremely 
high-velocity outflow created at the early post- 
AGB phase of stars, that is observed as water 
fountains. This may be due to gas ejecti on by 
binar y-s tar interactions iN ordhaus & Blackmail 



20061 ). iMen'shchikov et al. I (|2002( ) found that 



the massive dust torus of the Red Rectangle 
has M ~ 1.2 Mq, which was formed in the 
ejection of a common envelope after Roche lobe 
overflo w. Similarly, stellar merging (for exam- 



ple, iBallv fc Zinnecker 1 120051 : iTvlenda fc Soker 



20061 ) may create such massive circumstellar ma- 
terial. Alternatively these stars may simply be in 
a dense dust cloud, as found in IRAS 19312-f 1950 
(jPeguchi et al. 1 l2004al) . We cannot deny such a 
possibility for some of these sources (but not all 
of them) because they are seen toward thick dust 
clouds. We infer here that massive circumstellar 
matter of the thickest OH/IR stars in the present 
sample can be formed by one of the aforemen- 
tioned mechanisms. 
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Institute of Technology, funded by the National 
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4. Conclusion 

SiO maser and OH/IR objects with no 2MASS 
counterparts were identified in the GLIMPSE 
3.6/4.5/5.8/8.0 /im images with a spatial reso- 
lution of 1.2". Searches for SiO and II2O masers 
lead to 10 new detections, each of which show 
characteristic maser spectra of the stronger SiO 
J = 1-0 z) = 2 to w = 1 line or wide spread II2O 
maser emission. It turned out that the dust en- 
velopes of these objects are exceptionally thick, 
though they suffer from a large interstellar extinc- 
tion. This suggests that a considerable mass of 
materials toward the objects must exist in the en- 
velope. The SiO detection rate was uncorrelated 
with the mass of the envelope, but the lower de- 
tection rate was obtained for the lower-expansion 
velocity subsample of the OH/IR stars. The latter 
suggests that the higher-expansion velocity ob- 
jects are AGB stars and lower-expansion velocity 
objects are more evloved stars such as post-AGB 
stars. The derived envelope mass is not correlated 
with this sequence of stellar evolution. From these 
observations, we believe that the thick material in 
some of these objects cannot be ascribed to the 
AGB-phase wind of the central star. 



8 



A. IRAS 19027+0517 



The Glimpse images of the SiO maser source, IRAS 19027+0517, exhibit conspicuous nebulosity around 
the central star at 5.8 and 8.0 fj,m (Fig. 6). The nebula is lopsided to the central star and extended to the 
north. The surface brightness of the nebula seems to increase toward the central star, which suggests a 
physical association of this nebula to the central star. Among 200 SiO sources investigated, only this object 
has an extended nebulosity. The 2MASS images show a medium bright central star {K = 8.2), but do not 
exhibit any extended features. The nebulous feature is cons idered to be due to the PAH emission with sharp 
spectral features at 6.2 and 7.8 /xm ()Verstraete et al." 2001 ). It is quite unusual that the 0-rich central star 
in the AGB phase (with SiO masers) accompanies such extended PAH nebulosity. Therefore, we particularly 
investigated this object more with various molecular lines with the NRO 45m telescope. In addition to the 
previous detections of the SiO J = 1-0 v — 1 and 2 lines (jDeguchi et al. I l2004bh . we found emissions of the 
SiO J = 1-0 V = and 3, and J = 2-1 v = 1 lines; the spectra are shown in figure 7. Furthermore, we made 
45"-grid 5-point mapping observations by CO, HCN and HCO+ J = 1-0 lines, but no emission associated 
with the central star was detected, except for the fore/background clouds with large velocity separations 
(see figure 8). The observational results are summarized in table 5. Therefore, we conclude that most of 
the molecules are dissociated in the nebula toward IRAS 19027+0517, an 0-rich AGB star, by UV radiation 
from a closeby UV source creating the lopsided nebulosity. 
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30 


9 


-04 


34 


35 


7.1 





081 


-0.207 


n 


n — 





^"y" or "n" indicates the detection or nondetect ion i n this paper, o thewi se no ted. References : ^: [ Nalcash ima &: De^uchi 
^: lEngcls ot al. I 1119861). '^i ISevenster et al l 1 120011) . IGomoz et al.l 119901). ^: llzumiura et al.l Ill99gfl. " : " ' 

[Nyman ct al. 1 1998')'~ ^: ILikkel et al.l lll992f). lEngcTs & Lewis 1 1996^. Ite LinteLHekkert et al. I II1991I) . 
iLewis ct al.i U987i) . IChengalur et al.l jl993) . '^i iDeacon et al. i yOOj). 



112003b 
19861). 



lEder et alT II1988I) 
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Table 2 

Objects identified with Spitzer-Glimpse Survey. 



Source 


I 


6 


Fs.e F4..5 


F5.8 Fa.o [3.6] 


[3.6 


Ak 


fa 


(IRAS or MSX) 


(°) 


(°) 


(mJy) (mJy) 


(mJy) (mJy) -[5.8] 


-[12 







18034—1858 


+ 11.0061 


+0.9222 


29.9 


308.8 


1617.0 




5.30 


8.74 


1.50 


7.6 


18100 — 1915 


+ 11.5216 


— 0.5824 


72.7 


712.7 


3641.0 




5.21 


8.21 


3.05 


4.9 


18135 — 1456 


+ 15.7011 


+0.7706 


26.2 


166.6 


1033.0 




4.95 


10.06 


1.74 


15.9 


18161 — 1713 


+ 13.9884 


— 0.8611 


264.8 


1988.0 


11510.0 




5.06 


7.52 


3.33 


3.2 


18182 — 1504 


+ 16.1173 


—0.2907 


512.4 




14120.0 




4.57 


7.51 


5.64 


2.6 


18198—1249 


+18.2956 


+0.4291 


56.2 


847.5 


1818.0 




4.74 


8.20 


2.49 


5.1 


18212 — 1227 


+18.7688 


+0.3017 


1.9 


9.9 


77.7 


562.0 


4.99 


10.76 


3.20 


20.9 


18241-1010 


+21.1166 


+0.7775 


676.4 


2918.0 


9559.0 




3.84 


5.73 


1.65 


1.3 


18245-1248 


+18.8387 


-0.5622 


81.2 


447.1 


1468.0 


1750.0 


4.11 


5.40 


4.91 


0.8 


18257-1000 


+21.4565 


+0.4911 


434.1 










5.93 


2.63 


1.4 


18286-0959 


+21.7972 


-0.1272 


541.6 










7.37 


7.71 


2.0 


18298-0904 


+22.7483 


+0.0252 


231.0 




3409.0 




3.89 


7.31 


7.15 


2.1 


18325-0721 


+24.5814 


+0.2243 


198.1 


600.7 


2845.0 




3.86 


6.01 


11.48 


0.3 


18327-0645 


+25.1301 


+0.4842 


42.9 


270.2 


1074.0 




4.46 


8.50 


4.09 


5.3 


18407-0358 


+28.5203 


-0.0143 


75.3 




2288.0 




4.67 


6.82 


10.74 


1.1 


18424-0509 


+27.6622 


-0.9174 


98.3 


463.1 


1573.0 




3.98 


7.99 


1.31 


5.0 


18444-0359 


+28.9303 


-0.8288 


232.9 




3738.0 




3.98 


6.38 


1.77 


1.9 


18450-0148 


+30.9439 


+0.0351 


1.0 


12.5 


0217.2 


2564.0 


6.81 


13.51 


14.79 


37.8 


18460-0151 


+31.0124 


-0.2194 


23.9 


374.8 


2598.0 




6.06 


9.54 


8.23 


6.7 


18460-0254 


+30.0910 


-0.6865 


558.9 


5860.0 








8.47 


3.35 


5.5 


18488-0107 


+31.9845 


-0.4853 


778.6 


6270.0 


18690.0 




4.42 


6.57 


2.70 


2.0 


18509-0018 


+32.9528 


-0.5689 


273.0 


1198.0 


3568.0 


3952.0 


3.76 


7.47 


1.59 


3.6 


18517+0037 


+33.8727 


-0.3353 


107.9 




4579.0 




5.03 


8.52 


3.62 


5.6 


18596+0315 


+37.1184 


-0.8474 


9.9 


49.7 


217.3 


806.0 


4.32 


8.73 


1.54 


7.5 


19006+0624 


+40.0224 


+0.3813 


85.7 


484.3 


1566.0 




4.12 


6.09 


2.76 


1.5 


19011+0638 


+40.2794 


+0.3761 


84.9 


939.7 


5077.0 




5.41 


7.89 


2.71 


4.2 


19025+0702 


+40.8006 


+0.2452 


22.1 


187.6 


1009.0 


4079.0 


5.11 


8.80 


2.54 


7.2 


19087+1006 


+44.2406 


+0.3086 


00.9 


14.7 


146.5 


971.8 


6.49 


11.37 


3.80 


28.3 


19112+1220 


+46.4995 


+0.8093 


01.6 


23.5 


193.8 


1029.0 


6.17 


11.35 


1.88 


32.9 


19114+0920 


+43.8677 


-0.6248 


192.3 


865.6 


2582.0 


1809.0 


3.78 


6.91 


2.76 


2.4 


19180+1230 


+47.4258 


-0.5625 


93.5 


513.5 


1714.0 




4.12 


7.10 


2.92 


2.6 


19254+1631 


+51.8039 


-0.2249 


155.1 


1076.0 


4806.0 




4.69 


7.95 


2.80 


4.3 


19254+1724 


+52.5815 


+0.2014 


0.2 


02.4 


44.4 


503.2 


6.83 


13.14 


2.40 


88.4 


19266+1815 


+53.4611 


+0.3551 


42.8 


253.3 


910.0 


606.5 


4.28 


8.35 


2.03 


5.8 


19275+1702 


+52.5040 


-0.4085 


359.6 


1483.0 


4437.0 




3.69 


6.10 


2.85 


1.5 


19440+2251 


+59.4784 


-0.8966 


179.4 


1460.0 








7.53 


0.94 


4.0 


G014.2758 


+ 14.2760 


-0.1600 


367.1 










5.62 


8.35 


0.4 


G017.3913 


+ 17.3915 


-0.2892 


83.4 


583.4 


3113.0 




4.89 


7.15 


5.82 


2.1 


G022.0425 


+22.0428 


-0.6085 


6.6 


22.0 


83.5 


365.7 


3.72 


8.22 


3.91 


4.6 


G024.6610 


+24.6608 


+0.0867 


6.1 


130.0 


990.5 




6.49 


9.33 


11.05 


4.7 


G027.8433 


+27.8430 


-0.1061 


236.1 


811.8 


2035.0 


2720.0 


3.30 


6.26 


6.30 


1.2 



13 



Table 3: Detections by the SiO or H2O maser line. 



kjource 


Itansition 




(km s^^) 




Ta* 


line tlux 
km s^i") 




rms 
(K) 


i«inn— 1Q1 "i 


kJlW (J — X u c — 














.121 






2 


16 5 


\ 


.202 


9 087 


n 


.114 




I— 1 n f 1 (~\ 1 - — ^i-in 

n2\J U16 023 












n 

u . 


1 no 




kJlv7 f7 — -L u (y — 


n 


18 3 


n 




9 536 





1 53 




SiO r — 1-0 7) — 

oiv7 tJ — _L u (y — 


I 


22.9 


1 6 


.UU J- 


4^ 758 


n 


1 66 




SiO J = l-O V = 


2 


22^5 


41 


.500 


130.23 


0, 


.168 




SiO J = 1-0 = 


3 


22.5 


10 


.434 


28.839 


0, 


.173 






_ 1 
— ± 




1 
± 


.uyo 






,101. 


1 »1 Q8_1 24Q 


SiO 7 — 1-fl 7) — 


I 










n 


. J- Urt 




cliO 7^ — 1-0 71 — 


9 




n 

u 




n 49Q 





107 
,±\j i 


18241 — inin 


filO 7 — 1-0 ?? — 

kjivy U — ± \J u — 


\ 













• \Jijrt 




^iO 7^ — 1-0 7) — 


9 


J. J-U.O 











• VJOO 


18257— innn 


Q\r\ T _ in _ 

VJLKJ U J. U C 


1 










0, 


.097 




SiO J = 1-0 V = 


2 


114.9 





.781 


1.427 


0, 


.097 


18286-0959 


SiO J = 1-0 V = 


1 










0, 


.123 




SiO J = 1-0 V = 


2 










0, 


.118 




HoO 6ii^ 5oQ 




-14.7 


3 


.863 


166.3 





043 




SiO J = 2 1 w = 


1 










0, 


.065 




HCO+ J = 1 












0, 


.068 


18407-0358 


SiO J = 1 t- = 


1 


106.5 





.572 


1.457 


0, 


.108 




SiO J = 1- w = 


2 


107.3 





.921 


2.126 


0, 


.097 




H20 616-523 












0, 


.087 


18460-0151 


SiO J = 1-0 = 


1 










0, 


.100 




SiO J = 1-0 = 


2 










0, 


.084 




H2O 6i6~523 




117.0 


4, 


.881 


24.38 


0, 


.050 




SiO J = 2-1 = 


1 










0, 


.069 




HCO+ J = 1-0 












0, 


.084 


G017.3913 


SiO J = 1 Ov = 


1 










0, 


.119 




SiO J = 1-0 w = 


2 


28.7 


1 


.398 


4.106 


0, 


.117 


G024.6610 


SiO J = 1-0 u = 


1 










0, 


.103 




SiO J = 1-0 w = 


2 


57.7 


1 


.247 


3.662 


0, 


.101 



Table 4: Negative results for the SiO and H2O maser line search. 



Source 


rms 


rms 


rms'* 




(SiO 7; = 1) 


(SiO v = 2) 


(H2O) 




(K) 


(K) 


(K) 


18034-1858 


0.11b 


0.114 


0.096 


18161-1713 


0.119 


0.115 


0.107 


18212-1227 


0.096 


0.095 


0.086 


18298-0904 


0.097 


0.105 


0.075 


18325-0721 


0.075 


0.073 


0.087 


18327-0645 


0.099 


0.097 


0.088 


18424-0509 


0.357 


0.396 


0.088 


18444-0359 


0.126 


0.121 


0.113 


18596+0315 


0.117 


0.115 


0.086 


19011+0638 


0.113 


0.102 


0.084 


19025+0702 


0.111 


0.113 




19087+1006 


0.077 


0.075 


0.080 


19112+1220 


0.088 


0.084 


0.082 


19114+0920 


0.086 


0.086 


0.151 


19180+1230 


0.071 


0.073 




19254+1724 


0.094 


0.082 




19266+1815 


0.077 


0.071 




19275+1702 


0.094 


0.090 


0.093 


G014.2758 


0.109 


0.109 


0.099 


G022.0425 


0.075 


0.078 




G027.8433 


0.103 


0.099 


0.062 



"Dash indicates not observed. 
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Tublo o: Obhcrvrttiuii suimiirtry lur IRAS li)()27+()517. 



I'ransition 




rest Ireq. 


Vlsr 


^ a 


lino Hux 


rms 






(GHz) (k 


m s~^) 


(K) 


(km K) 


(K) 


H2U 616 523 




22.28508 








0.078 


NH3 ll-ll 




28.694 








0.068 


NH3 22-22 




28.722 








0.044 


NH3 83-83 




28.870 








0.052 


SiO J = 1 Ov = 





48.42886 


84.2 


0.178 


0.417 


0.058 


SiO J = 10 V = 


1 


48.12208 


81.7 


1.288 


8.848 


0.058 


SiO J = 1-0 w 


2 


42.82059 


81.7 


1.202 


8.826 


0.067 


SiO J=l-Ov = 


8 


42.51938 


31.6 


0.298 


0.557 


0.053 


29SiO J = l-Ov 


= 1 


42.87992 








0.054 


H53q! 




42.95197 








0.062 


SiO J = 2-lv = 


1 


86.24342 


32.3 


0.204 


0.584 


0.034 


SiO J = 2-lv = 





86.84700 








0.037 


HCN J = 1-0 




88.63185 


a 






0.053 


HCO+ ,7=10 




89.18852 








0.061 


CO J = 1-0 




115.27120 


a 






0.200 


'Contaminations at VJ, 


sr 


12, 46, and 63 km s 


-1 
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1B444-O350 
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Fig. 1.— a. Composite-color images (285" x 285") of the GLIMPSE field of the som-ces in Table 1. IRAC 3.6, 
4.5 and 5.8 fim bands are represented by blue, green, and red colors respectively. Objects are at the center 
of the images, and the directions of increasing Galactic longitude and latitude are left and up, respectively. 
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Fig. 1. — b. and c. — Continued. 
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Fig. 2.— SiO J = 1-0 = 1 and 2 spectra of the detected sources. For IRAS 18182-1504 the SiO J = 1-0 
u = 1, 2, and 3, and ^^SiO J = 1-0 v = Q spectra are shown. Source name and observed date (in yymmdd.d 
format) are indicated on the left of each panel, were indicated 
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Fig. 3. — H2O 6i6-523 spectra for the detected sources. Source name and observed date (in yymmdd.d 
format) are indicated on the left of each panel. 
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Fig. 4. — Two-color diagram, [3.6] — [12] versus [3.6] — [5.8]. The filled and unfilled circles indicate SiO 
detection and nondetection, and the square indicates the position of OH 127.8+0.0. The dotted arrow 
indicates direction of interstellar reddening; the length corresponds to Ak = 15. 




Fig. 5. — Excess factor versus expansion velocity of the OH 1612 MHz double-peak objects. Filled and 
unfilled circles indicate SiO detection and nondetection. The broken line indicates a separation for the 
low-mass stars. 



Fig. 6.— Glimpse color-composite image (285" x 285") of IRAS 19027+0517 made from 4.5 (blue), 5.8(green), 
and 8.0 (red) /im IRAS image. The directions of increasing Galactic longitude and latitude are left and up, 
respectively. 
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Fig. 7.— SiO, H2O and NH3 spectra (a) CO J = 1-0 spectra (b) toward IRAS 19027+0517. The number 

between the parenthesis on the left indicates R.A. and Dec. offsets in arcsec from the star position. The top 
is the spectrum of the bottom (center) subtracted by the average of surrounding 4 positions (middle). 
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Fig. 8. — CO J = 1-0 spectra toward IRAS 19027+0517. The number between the parenthesis on the left 
indicates R.A. and Dec. offsets in arcsec from the star position. The top is the spectrum of the bottom 
(center) subtracted by the average of surrounding 4 positions (middle). 
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